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Large-scale  molecular  dynamics  (MD)  simulations  are  used  to  understand  the  macroscopic  yield 
behavior  of  nanocrystalline  Cu  with  an  average  grain  size  of  6  nm  at  high  strain  rates.  The  MD 
simulations  at  strain  rates  varying  from  10^  s“^  to  8  x  10^  s“^  suggest  an  asymmetry  in  the  flow 
stress  values  in  tension  and  compression,  with  the  nanocrystalline  metal  being  stronger  in 
compression  than  in  tension.  The  tension-compression  strength  asymmetry  is  very  small  at 
10^  s“\  but  increases  with  increasing  strain  rate.  The  calculated  yield  stresses  and  flow  stresses 
under  combined  biaxial  loading  conditions  (X-  Y)  gives  a  locus  of  points  that  can  be  described 
with  a  traditional  ellipse.  An  asymmetry  parameter  is  introduced  that  allows  for  the  incorpo¬ 
ration  of  the  small  tension-compression  asymmetry.  The  biaxial  yield  surface  {X-Y)  is  calculated 
for  different  values  of  stress  in  the  Z  direction,  the  superposition  of  which  gives  a  full  three- 
dimensional  (3-D)  yield  surface.  The  3-D  yield  surface  shows  a  cylinder  that  is  symmetric 
around  the  hydrostatic  axis.  These  results  suggest  that  a  von  Mises-type  yield  criterion  can  be 
used  to  understand  the  macroscopic  deformation  behavior  of  nanocrystalline  Cu  with  a  grain 
size  in  the  inverse  Hall-Petch  regime  at  high  strain  rates. 
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I.  INTRODUCTION 

Nanocrystalline  metals,  due  to  their 
enhanced  strength  and  wear  resistance,  show  significant 
promise  for  a  number  of  future  high-technology  appli¬ 
cations  such  as  microelectromechanical  systems,  nano¬ 
electromechanical  systems,  and  nanoscale  devices.  To 
maximize  the  performance  and  reliability  of  these 
devices,  a  fine-tuning  of  the  nanostructure  is  required, 
which  necessitates  a  fundamental  understanding  of  the 
deformation  and  failure  mechanisms. Nanocrystalline 
metals  have  been  extensively  studied, to  understand 
the  micromechanisms  governing  their  macroscopic 
mechanical  behavior.  At  the  atomic  scale,  plastic  defor¬ 
mation  mechanisms  in  nanocrystalline  materials  can  be 
classified  into  dislocation-based  and  grain-boundary 
(GB)-based  processes.  A  reduction  in  grain  size  results 
in  an  increase  in  the  yield  strength  of  materials,  a 
relation  known  as  the  Hall-Petch  effect, which  states 
that  the  yield  strength  varies  inversely  with  the  square 
root  of  the  grain  size.  Recent  studies^^^“^"^^  indicate  that 
the  increase  in  strength  with  decreasing  grain  size 
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reaches  a  maximum  after  which  a  further  decrease  in 
the  grain  size  (less  than  ~  12  nm)  can  result  in  weaken¬ 
ing  of  the  metal.  This  weakening  of  the  metal,  due  to  the 
shift  in  the  dominating  mechanism  of  plastic  deforma¬ 
tion  from  dislocation-induced  plasticity  in  the  case  of 
coarse-grained  materials  to  GB  sliding  in  the  case 
of  ultrasmall  grain  sizes,  is  referred  to  as  the  inverse 
Hall-Petch  behavior. 

Nanocrystalline  metals  with  ultrasmall  grain  sizes 
{d  <  30  nm)  have  gained  considerable  attention  due  to 
their  increased  strengths  during  deformation  at  high 
strain  rates. In  addition,  shock  loading  of  these 
ultrasmall  nanocrystalline  metals  at  speeds  that  are 
greater  than  the  speed  of  sound  (strain  rates  >10^  s“^) 
limits  the  GB  sliding  mechanism  and  thus  results  in 
ultrahigh  strength  values  for  the  nanocrystalline  me¬ 
tal.  This  limiting  of  GB  sliding  at  ultrahigh  strain 
rates  may  lead  to  modifications  in  the  yield  criteria/ 
behavior  at  the  ultrasmall  grain  sizes.  To  maximize  the 
performance  and  reliability  of  these  devices,  a  fine- 
tuning  of  the  nanostructure  is  required,  which  necessi¬ 
tates  a  fundamental  understanding  of  the  deformation 
and  failure  of  the  constituent  materials  of  the  devices. 

The  macroscopic  response  of  materials  to  plastic 
deformation  can  be  understood  by  calculating  the  yield 
criterion  governing  the  multiaxial  plasticity  at  high 
strain  rates.  The  most  commonly  used  phenomenolog¬ 
ical  yield  criteria^^^’^^^  for  polycrystalline  metals,  such  as 
von  Mises  and  Tresca,  are  based  on  the  assumption  that 
the  yielding  of  the  metal  is  determined  by  the  deviatoric 
stress  and  is  independent  of  the  pressure  and  the  normal 
stress.  These  conventional  yield  criteria  have  been  found 
to  be  appropriate  to  study  deformation  behavior  in 
polycrystalline  metals. Plasticity  in  polycrystalline 
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metals  is  dominated  by  the  nucleation  and  propagation 
of  dislocations.  Plasticity  in  nanocrystalline  metals, 
however,  has  contributions  from  dislocation-based  and 
GB-based  processes.  The  GB  processes  become  more 
dominant  at  grain  sizes  less  than  30  nm.  Particularly 
important  is  the  case  of  nanocrystalline  metals  with  the 
grain  size  in  the  inverse  Hall-Petch  regime  {d<\5  nm) 
in  which  GB  sliding  dominates  the  deformation  behav¬ 
ior. Deformation  studies  of  nanocrystalline  metals 
(2  to  4  nm)  using  molecular  static  simulations  indicate  a 
strong  asymmetry  in  tensile  and  compressive 
strengths. ^  A  similar  behavior  is  observed  in  metallic 
glasses,  in  which  a  strong  asymmetry  in  tension  and 
compression  is  observed  and  can  be  well  described  using 
the  Mohr-Coulomb  criterion,  which  accounts  for  a 
normal  stress  dependence. This  normal  stress- 
dependent  criterion  has  been  recently  extended  to  the 
ultrasmall  grains  (2  nm),  suggesting  that  ultrasmall 
grains  can  be  considered  as  amorphous  materials  at 
the  lower  limit  of  the  grain  size.^^"^^ 

To  date,  most  of  the  crystal  plasticity  models  devel¬ 
oped  to  describe  plastic  deformation  and  failure  in 
conventional  coarse-grained  metals  do  not  explicitly 
include  GB  effects  on  plasticity.  As  a  result,  these 
models  may  not  be  suitable  for  predicting  the  material 
response  at  grain  sizes  in  the  inverse  Hall-Petch  regime. 
The  modification  of  conventional  continuum  failure 
models  to  accurately  predict  material  behavior  at  the 
nanoscale  would  represent  a  paradigm  shift  in  the  design 
and  realization  of  revolutionary  nanocrystalline  materi¬ 
als.  The  molecular  dynamics  (MD)  simulation  technique 
has  the  potential  to  provide  the  atomic-level  structural 
information  about  the  deformation  mechanisms  needed 
to  modify  and  expand  the  capabilities  of  continuum 
simulations.  The  focus  of  this  article  is  to  understand  the 
response  of  these  ultrasmall-grain-size  nanocrystalline 
metals  during  deformation  under  conditions  of  multi- 
axial  loading.  The  strain  rates  used  here  are  chosen  to 
reproduce  those  observed  during  shock  loading  to  create 
ultrahigh  strengths  in  nanocrystalline  metals.^^^^  Shock 
loading  using  short  (nanosecond)  laser  pulses  leads  to 
peak  strain  rates  exceeding  10^  s“V^^  ^  whereas  the 

use  of  ultrashort  (femtosecond)  laser  pulses  results  in 
strain  rates  exceeding  10^  ^-i  [28,29]  ^  result,  MD 

simulations  are  carried  out  to  understand  the  macro¬ 
scopic  deformation  behavior  of  nanocrystalline  Cu  with 
an  average  grain  size  of  6  nm  at  ultrahigh  strain  rates 
(>10^  s“^).  Three  aspects  of  deformation  behavior  are 
studied:  the  tension-compression  strength  asymmetry, 
the  biaxial  yield  surface,  and  the  three-dimensional 
(3-D)  yield  surface.  The  computational  details  are 
presented  in  Section  II.  The  tension-compression 
strength  asymmetry  is  discussed  in  Section  III,  followed 
by  the  calculation  of  the  3-D  yield  surface  in  Section  IV. 


11.  COMPUTATIONAL  METHODS 

The  initial  MD  system  of  nanocrystalline  Cu  used  in 
this  study  is  created  by  the  Voronoi  construction 
method,  as  suggested  by  Derlet  and  Van  Swygenho- 
ven.^^^^  Periodic  boundary  conditions  are  used  in  all 


three  directions.  In  addition,  the  system  is  equilibrated 
at  300  K  for  200  picoseconds  to  relax  the  structure  after 
the  Voronoi  construction.  The  density  of  the  final 
configuration  is  calculated  to  be  ~99  pet  of  the  bulk 
density  of  copper.  This  procedure  constructs  a  system 
with  random  grain  orientations  containing  no  textures 
and  a  grain  size  distribution  close  to  a  log-normal 
distribution.^^  The  centro  symmetry  parameter 
(Aesp)^^^^  is  calculated  for  each  atom  as 

Aesp  =  ^  1^/  +  [1] 

/=1,6 

The  blue  color  (Aesp  ~  0)  corresponds  to  FCC  stack¬ 
ing;  the  green  color  (Aesp  ~  8)  represents  grain  bound¬ 
aries  and  partial  dislocations;  and  the  red  color 
(Aesp  >15)  represents  a  surface.  A  typical  system  with 
an  average  grain  size  of  6  nm  and  containing  122 
grains  is  shown  in  Figure  1,  with  the  atoms  colored 
according  to  the  Aesp  values.  The  system  consists  of 
approximately  1.2  million  atoms  arranged  in  122  grains 
with  an  average  grain  diameter  of  ~6  nm.  The  MD 
simulations  are  carried  out  using  the  using  the  Voter- 
Chen  (VC)  formulation^^^^  of  the  embedded  atom 
method  potential  for  copper.  The  VC  potential  pro¬ 
vides  a  good  description  of  the  unstable  and  stable 
stacking  fault  energies  and  therefore  is  well  suited  to 
describe  deformation  behavior  for  Cu.^^"^^ 


0  2  4  6  8  10  12  14  16 


Fig.  1 — Initial  configuration  of  nanocrystalline  Cu  with  an  average 
grain  size  of  6  nm.  System  consists  of  approximately  1.2  million 
atoms  arranged  in  122  grains;  each  atom  is  colored  according  to  the 
^esp- 
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The  elements  of  the  tensor  of  atomic-level  stresses  are 
calculated  as 


L2o 


9  IJ  IJ  ^  ^  I  I 


[2] 


where  a  and  j]  are  the  Cartesian  components,  Qq  is  the 
atomic  volume,  Fg  is  the  force  on  atom  i  due  to  atom 
j\  Mi  is  the  mass  of  atom  /,  and  is  the  velocity  of 
atom  i.  An  effective  stress  (a^)  and  an  effective  strain 
according  to  the  von  Mises  criterion  for  yielding  is 
calculated  as 


(7^ 


6  at 


-)) 


1  1/2 


[3] 


[4] 


where  Gy,  and  g^  are  the  stresses  averaged  over  the 
entire  system  in  the  X,  Y,  and  Z  directions,  respectively, 
and  Sy,  and  are  the  engineering  strains  in  the  X,  Y, 
and  Z  directions,  respectively.  The  time-step  for  all  of 
the  MD  simulation  runs  is  chosen  to  be  0.002  ps.  The 
temperature  is  allowed  to  evolve  during  the  deformation 
process.  The  effect  of  ultrahigh  strain  rate  deformation 
on  the  strength  asymmetry  is  discussed  in  Section  III 
along  with  the  predicted  plastic-yield  surface. 


III.  TENSION-COMPRESSION  STRENGTH 
ASYMMETRY  IN  NANO  CRYSTALLINE  Cu 

The  MD  simulations  of  uniaxial  loading  of  the  nano¬ 
crystalline  Cu  sample  are  carried  out  at  constant  strain  by 
deforming  the  sample  in  the  X  direction,  and  maintaining 
zero  stress  conditions  in  the  Y  and  Z  directions 
(cr^  =  (jz  =  0  and  g^  =  g).  At  each  time-step,  deforma¬ 
tion  is  achieved  by  adjusting  the  x  coordinate  of  all  the 
atoms  and,  accordingly,  the  length  of  the  system,  in  the  X 
direction  by  a  scaling  parameter  chosen  to  achieve  the 
desired  strain  rate  (e).  The  coordinates  in  the  Y  and  Z 
directions  are  scaled  so  as  to  have  zero  stress  in  the  Y  and 
Z  directions.  The  MD  simulations  of  deformation  are 
carried  out  for  strain  rates  ranging  from  1  x  10^  s“^  to 
8  X  10^  s“^  The  choice  of  strain  rates  is  based  on  those 
generated  during  shock  loading  using  picosecond  and 
femtosecond  laser  pulses. Based  on  Eqs.  [3]  and 
[4],  for  the  conditions  of  uniaxial  stress  loading,  the 
effective  von  Mises  stress  (g^)  reduces  to  the  stress  in  the 
loading  direction  (g^),  and  the  effective  strain  reduces 
to  the  strain  in  the  loading  direction  The  plots  of  the 
effective  von  Mises  stress  as  a  function  of  effective  strain 
during  uniaxial  deformation  in  tension  in  comparison 
with  that  in  compression  are  shown  for  strain  rates  of 
1  X  lO**  s“'  to  8  X  10^  s  ^  in  Figure  2.  The  curves  are 
initially  linear  up  to  the  yield  point,  after  which  they  start 


Fig.  2 — Plots  of  effective  von  Mises  stress  as  a  function  of  effective 
strain  in  tension  and  compression  at  strain  rates  of  1  x  10^  s“^  and 
8  X  10^  s“\  respectively.  Nanocrystalline  Cu  shows  greater  strength 
in  compression  (C)  as  compared  to  tension  (T). 


Table  I.  Values  of  Flow  (Peak)  Stress  during  Uniaxial 
Deformation  of  Nanocrystalline  Cu  with  an  Average  Grain 
Size  of  6  Nanometers  at  Various  Strain  Rates  in  Tension 
and  Compression 


Strain 

Rate  s  (s“^) 

Tension 
(Flow  Stress) 
ffj  (GPa) 

Compression 
(Flow  Stress) 
(GPa) 

Strength 

Asymmetry  (Pet) 

1  X  10® 

1.986 

2.101 

5.79 

2  X  10® 

2.164 

2.303 

6.42 

4  X  10® 

2.405 

2.577 

7.15 

8  X  10® 

2.794 

3.010 

7.73 

1  X  10" 

6.399 

8.712 

36.15 

to  deviate  from  elastic  behavior.  The  flow  stress  (oA  is 
defined  as  the  peak  stress  in  the  stress-strain  curve. 

The  nanocrystalline  Cu  is  observed  to  be  stronger  in 
compression  than  in  tension.  This  strength  asymmetry, 
defined  as  the  tension-compression  strength  differen¬ 
tial,  is  larger  at  the  higher  strain  rates.  The  calculated 

values  of  the  flow  stress  in  tension  and  compression  at 
various  strain  rates  are  tabulated  in  Table  I  and  plotted  in 
Figure  3.  The  nanocrystalline  sample  shows  a  strength 
asymmetry  of  ~8  pet  at  8  x  10^  s“\  which  is  much  lower 
than  that  (~35  pet)  observed  in  nanocrystalline  Ni  during 
quasi-static  loading  using  molecular  static  simulations. 

It  should  be  noted  that  the  quasi-static  deformation  is 
obtained  by  scaling  the  atomic  coordinates  (strain  incre¬ 
ments  along  the  axis)  followed  by  relaxation  of  the  system 
using  the  conjugant  gradient  technique. This  relax¬ 
ation  in  molecular  static  simulations,  however,  does  not 
allow  for  complete  adiabatic  relaxation  of  the  system 
during  the  deformation  process.  The  deformation  strains 
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Fig.  3 — Plot  of  the  calculated  yield  stress  and  flow  (peak)  stress  in 
tension  and  compression  at  strain  rates  1  x  10^  s“^  to  8  x  10^  s“\ 
respectively.  Yield  stress  and  flow  stress  values  are  greater  in  com¬ 
pression  and  the  difference  increases  with  increasing  strain  rates. 


used  may  therefore  result  in  deformation  at  ultrahigh 
strain  rates,  which  results  in  the  large  strength  asymmetry 
(~35  pet)  in  tension  and  compression  for  nanocrystalline 
Ni.  To  test  this,  the  flow  stress  is  calculated  for  uniaxial 
stress  loading  in  tension  and  compression  using  MD 
simulations  at  a  strain  rate  of  10^^  s“^  The  tension- 
compression  strength  asymmetry  is  calculated  to  be 
~36  pet  at  a  strain  rate  of  10^^  s“^  This  value  of  the 
asymmetry  compares  very  well  with  that  observed  for 
quasi-static  loading  of  nanocrystalline  Ni  and  suggests 
that  the  large  asymmetry  observed  in  quasi-static  defor¬ 
mation  of  Ni  using  molecular  static  simulations  may  be 
due  to  the  unusually  high  strain  rates  of  deformation 
induced  by  the  technique. 


IV.  PLASTIC-YIELD  SURFACE 

A  full  3-D  plastic  yield  surface  using  MD  simulations 
can  be  obtained  by  deforming  the  nanocrystalline  metal 
by  equal/unequal  amounts  in  the  X,  Y,  and  Z  directions 
and  plotting  the  plastic  yield/flow  stresses  in  the  X,  Y, 
and  Z  axes.  While  the  process  can  be  complicated  and  in 
turn  tedious,  the  3-D  yield  surface  can  also  obtained  by 
computing  the  biaxial  yield  surfaces.  Here,  in  this  work, 
the  biaxial  yield  surface  is  be  calculated  by  plotting  the 
yield/flow  stresses  during  loading  of  the  nanocrystalline 
metal  by  equal/unequal  amounts  in  the  X  and  Y 
directions  and  keeping  the  stress  in  the  Z  direction 
constant  =  0  GPa).  The  biaxial  yield  surface  is  be 
then  calculated  in  a  similar  way  for  different  values  of 
the  constant  stress  in  the  Z  direction  (cr^  =  1  GPa; 
(Tz  =  —I  GPa).  The  calculated  biaxial  yield  surfaces  are 
then  be  superimposed  along  the  Z  axis  to  obtain  a  3-D 
yield  surface. 


Fig.  4 — Evolution  of  the  stresses  in  the  X,  Y,  and  Z  directions,  along 
with  the  effective  von  Mises  stress  as  a  function  of  effective  strain  for 
unbalanced  biaxial  loading  (X  in  tension  and  Y  in  compression 
shown).  Yield  stresses  in  X  and  Y  directions  are  calculated  based  on 
the  0.2  pet  offset  stress  on  the  von  Mises  curve.  Flow  stress  in  the  X 
and  Y  directions  is  calculated  at  the  value  of  the  effective  strain  cor¬ 
responding  to  a  peak  in  the  von  Mises  curve.  Effective  strains  corre¬ 
sponding  to  the  yield  stress  and  flow  stress  are  shown  by  the  dashed- 
dotted  lines. 

A.  Biaxial  Yield  Surface  at  =  0  GPa 

The  initial  nanocrystalline  Cu  system  is  equilibrated 
at  zero  pressure  for  50  ps  to  have  a  stress  of  cr^  =  0  GPa 
in  the  Z,  T,  and  Z  directions.  To  calculate  the  biaxial 
yield  surface,  deformation  is  achieved  at  each  time-step 
by  adjusting  the  x  and  y  coordinates  of  all  the  atoms  by 
a  scaling  parameter  based  on  the  loading  conditions  in 
the  X  and  Y  directions.  The  scaling  in  the  X  and  Y 
directions  may  be  equal  or  unequal  depending  on  the 
desired  combination  of  stresses  desired  =  Gi, 
(Ty  =  (72,  and  Gz  =  G3  =  0).  The  maximum  scaling  used 
gives  a  strain  rate  of  10^  s“^  For  the  case  of  unbalanced 
deformation  in  the  X  and  Y  directions,  the  calculation  of 
yield  stresses  and  flow  stresses  can  be  complicated.  For 
example.  Figure  4  shows  the  evolution  of  the  stresses 
averaged  in  the  X,  Y,  and  Z  directions,  along  with  the 
von  Mises  stress  as  a  function  of  effective  strain.  The 
deformation  conditions  correspond  to  tensile  loading  in 
the  X  direction  and  compressive  loading  in  Y  direction. 
Yielding  is  defined  by  the  0.2  pet  offset  stress  in  each 
direction,  as  shown  by  the  black  lines  in  Figure  4.  It  can 
be  seen  that  the  yielding  on  the  effective  stress  curve 
occurs  at  the  same  value  of  effective  strain  as  observed 
for  yielding  in  the  X  and  Y  directions.  As  a  result,  for  the 
various  combinations  of  loading  conditions,  the  yield 
stress  values  in  the  X  direction  (cr^)  and  Y  direction  (o^) 
are  calculated  at  the  effective  strain  corresponding  to  the 
0.2  pet  offset  stress  on  the  effective  stress  curve  (a^). 
Similarly,  the  flow  stress  (peak  value  of  the  stress)  values 

in  the  X  direction  [g{^  and  Y  direction  are 

calculated  at  the  effective  strain  corresponding  to  the 
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peak  value  of  the  effective  stress  curve  [gQ  .  The  effective 
strains  corresponding  to  the  yield  stress  and  flow  stress 
are  shown  by  the  dashed-dotted  lines  in  Figure  4. 
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Fig.  5 — Plot  of  the  locus  of  {a)  the  yield  stress  (0.2  pet  offset  stress) 
values  and  {b)  flow  stress  (peak  value  of  stress)  values,  for  biaxial 
loading  in  the  X  and  Y  directions  and  zero  stress  in  the  Z  direction. 
Dotted  ellipse  fitted  to  the  data  suggests  a  von  Mises-type  yield  sur¬ 
face.  Solid  lines  indicate  the  shifted  center  of  the  ellipse  to  account 
for  the  tension-compression  strength  asymmetry. 


The  locus  of  the  yield  stress  values  and  flow  stress 
values  for  biaxial  loading  =  Oi,  Oy  =  a 2,  and 
(73  =  =  0)  is  shown  in  Figures  5(a)  and  (b),  respec¬ 

tively.  The  yield  and  flow  surfaces  can  be  fitted  to  an 
ellipse  (dotted),  as  shown,  and  suggests  a  von  Mises-type 
surface  for  yielding.  The  center  of  the  ellipse,  however,  is 
shifted  by  a  very  small  amount  in  the  third  quadrant  in 
order  to  fit  the  yield  surface.  This  shift  in  the  center  of 
the  ellipse  is  due  to  the  inherent  tension-compression 
strength  asymmetry  at  these  high  strain  rates.  The 
amount  of  shift  in  the  center  of  the  ellipse  is  defined  as 
an  asymmetry  parameter  (a)  and  is  symmetric  in  the  X 
and  Y  directions.  The  traditional  von  Mises  ellipse  has 
the  minor  and  major  axis  lengths  defined  as 
and  (^Vl^Gy,  respectively,  where  Gy  is  the  calculated  yield 
stress  in  tension.  The  yield  stress  and  flow  stress  values 
for  the  fitted  ellipse  are  calculated  from  the  shifted 
center  of  the  ellipse,  as  shown  in  Figures  5(a)  and  (b), 
respectively.  The  parameters  of  the  fitted  ellipse  for  the 
yield  stresses  give  the  length  of  the  minor  axis  to  be 
1.05  GPa,  which  compares  very  well  with  the  value  of 
-1.06  GPa  predicted  by  a  traditional  von  Mises  ellipse 
using  Gy.  The  length  for  the  major  axis  for  the  fitted 
ellipse  is  calculated  to  be  2.03  GPa,  which  is  slightly 
larger  than  the  value  of  -1.84  GPa  predicted  by  a 
traditional  von  Mises  ellipse  using  Gy.  Similarly,  the  flow 
stress  ellipse  has  a  minor  axis  of  1.64  GPa,  which  agrees 
very  well  with  the  value  of  -1.67  GPa  predicted  by  a 
traditional  von  Mises  ellipse  using  Gf.  The  length  of  the 
major  axis  for  the  flow  stress  ellipse  is  calculated  to  be 
3.19  GPa,  which  is  slightly  larger  than  the  value  of 
-2.9  GPa  predicted  by  a  traditional  von  Mises  ellipse 
using  Gf.  Thus,  the  ellipse  fitted  to  the  yield  locus 
suggests  the  validity  of  the  von  Mises-type  yield  crite¬ 
rion,  but  deviates  slightly  in  the  length  of  the  major  axis 
or,  in  other  words,  the  aspect  ratio  of  the  predicted 
ellipse.  A  similar  fitting  can  be  done  for  the  Tresca 
criterion.  The  calculated  lengths  of  the  major  and  minor 
axis  are  tabulated  in  Table  II,  along  with  the  asymmetry 
parameter. 


B.  Biaxial  Yield  Surface  at  g^  =  1  GPa 

To  calculate  the  yield  stresses  at  a  stress  of  1  GPa  in 
the  Z  direction,  the  system  is  first  equilibrated  for  100  ps 
at  a  pressure  of  P  =  —1  GPa.  As  a  result,  the  stresses  in 
the  X,  F,  and  Z  directions  are  equilibrated  at  1  GPa 
prior  to  deformation.  The  deformation  is  then  carried 
out  by  scaling  the  x  and  y  coordinates  of  the  atoms  by 
equal  or  unequal  amounts,  depending  on  the  desired 
combination  of  loading  conditions  in  the  X  and  Y 


Table  II.  Calculated  Lengths  of  Major  and  Minor  Axes  of  Ellipse  Used  to  Fit  Biaxial  Yield  and  Flow  Surfaces  at  Various 
Constant  Values  of  Stress;  Gy  and  are  Yield  Stress  and  Flow  Stress  Values,  Respectively,  in  Tension  with  Shifted  Center 
of  Ellipse;  Xy  and  a/ Are  the  Shift  in  the  Center  of  the  Ellipse  from  the  Ideal  Value 


Z  Stress  (GPa) 

Yield  Stress 

Flow  Stress 

Minor  Axis 

Major  Axis 

Minor  Axis 

Major  Axis 

a/ 

Gz  =  0 

1.05 

2.03 

1.30 

0.04 

1.64 

3.19 

2.05 

0.06 

Gz  =  +1 

1.075 

2.15 

1.38 

0.12 

1.60 

3.10 

2.00 

0.05 

Gz  =  —1 

1.2 

2.25 

1.50 

0.10 

1.66 

3.35 

2.09 

0.06 
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Fig.  6 — Evolution  of  the  stresses  in  the  X,  Y,  and  Z  directions,  along 
with  the  effective  von  Mises  stress  as  a  function  of  the  effective 
strain  for  unbalanced  biaxial  loading  (X  in  compression  and  Y  in 
tension  shown).  Yield  stresses  in  X  and  Y  directions  are  calculated 
based  on  the  0.2  pet  offset  stress  on  the  von  Mises  curve.  Flow  stress 
in  the  X  and  Y  directions  is  calculated  at  the  value  of  effective  strain 
corresponding  to  a  peak  in  the  von  Mises  curve.  Effective  strains 
corresponding  to  the  yield  stress  and  flow  stress  are  shown  by  the 
dashed-dotted  lines. 


directions  =  di,  Gy  =  (72).  The  z  coordinates  of  the 
atoms  are  scaled  so  as  to  maintain  a  constant  stress  of 
(73  =  =  1  GPa  in  the  Z  direction.  Figure  6  illustrates 

the  offset  and  flow  stress  definitions  for  unbalanced 
deformation  with  compression  in  X  and  tension  in  Y.  It 
can  be  seen  that  the  stresses  in  the  X  and  Y  directions  are 
initially  at  1  GPa  and  evolve  linearly  during  the  defor¬ 
mation  prior  to  yielding.  The  stress  in  the  Z  direction  is 
held  constant  at  1  GPa.  The  calculated  values  of  the 
yield  stresses  and  the  flow  stresses  in  the  X  and  Y 
directions  for  various  loading  combinations  are  plotted 
in  Figures  7(a),  and  (b),  respectively.  The  yield  surface  is 
similar  to  that  calculated  at  zero  pressure,  and  can  be 
well  described  using  the  von  Mises-type  ellipse.  The 
center  of  the  ellipse  is  shifted  from  the  ideal  value  (1,1) 
by  a  very  small  amount  a  in  the  third  quadrant  due  to 
the  inherent  asymmetry  in  tension  and  compression.  The 
yield  stress  {oy)  and  flow  stress  (oy)  values  for  the  fitted 
ellipse  are  calculated  from  the  shifted  center  of  the 
ellipse,  as  shown  in  Figures  7(a)  and  (b),  respectively. 
The  lengths  of  the  major  and  minor  axes  for  the  fitted 
ellipse  are  tabulated  in  Table  II  along  with  the  shift  in 
the  center  of  the  ellipse  for  the  yield  and  flow  surfaces. 
The  fitting  parameters  give  the  length  of  the  minor  axis 
to  be  1.075  GPa,  which  compares  very  well  with  the 
value  of  -1.13  GPa  predicted  by  a  traditional  von  Mises 
ellipse  using  Oy.  Similarly,  the  length  for  the  major  axis 
for  the  yield  stress  ellipse  is  calculated  to  be  2.15  GPa, 
which  is  slightly  larger  than  the  value  of  -1.95  GPa 
predicted  by  a  traditional  von  Mises  ellipse  using  Oy. 
Similarly,  the  flow  stress  ellipse  has  a  minor  axis  of 
1.60  GPa,  which  agrees  very  well  with  the  value  of 


o; ,  (GPa) 

(a) 


(b) 

Fig.  7 — Plot  of  the  locus  of  {a)  yield  stress  (0.2  pet  offset  stress) 
values  and  {b)  flow  stress  (peak  value  of  stress)  values  for  biaxial 
loading  in  the  X  and  Y  directions  and  constant  stress  o^  0^=  \  GPa. 
Dotted  ellipse  fltted  to  the  data  suggests  a  von  Mises-type  yield  sur¬ 
face.  Solid  lines  indicate  the  shifted  center  of  the  ellipse  to  account 
for  the  tension-compression  strength  asymmetry. 


-1.63  GPa  predicted  by  a  traditional  von  Mises  ellipse 
using  Gf.  The  length  of  the  major  axis  for  the  flow  stress 
ellipse  is  calculated  to  be  3.10  GPa,  which  is  slightly 
larger  than  the  value  of  -2.83  GPa  predicted  by  a 
traditional  von  Mises  ellipse  using  oy. 

C.  Biaxial  Yield  Surface  at  g^  =  —1  GPa 

To  calculate  the  yield  stresses  at  a  stress  of  -1  GPa  in 
the  Z  direction,  the  system  is  first  equilibrated  for  100  ps 
at  a  pressure  of  P  =  1  GPa.  As  a  result,  the  stresses  in 
the  X,  F,  and  Z  directions  are  equilibrated  at  -1  GPa 
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Fig.  8 — Plot  of  the  locus  of  {a)  yield  stress  (0.2  pet  offset  stress)  val¬ 
ues  and  {b)  flow  stress  (peak  value  of  stress)  values  for  biaxial  load¬ 
ing  in  X  and  Y  directions  and  constant  stress  of  =  —  1  GPa. 
Dotted  ellipse  fitted  to  the  data  suggests  a  von  Mises-type  yield  sur¬ 
face.  Solid  lines  indicate  the  shifted  center  of  the  ellipse  to  account 
for  the  tension-compression  strength  asymmetry. 


prior  to  deformation.  Deformation  is  then  carried  out  as 
discussed  earlier.  The  locus  of  the  yield  and  flow  stresses 
in  the  X  and  Y  directions  at  =  - 1  GPa  is  shown  in 
Figures  8(a)  and  (b),  respectively.  The  calculated  values 
of  the  yield  stresses  and  flow  stresses  in  the  X  and  Y 
directions  for  various  loading  combinations  are  plotted 
in  Figures  8(a),  and  (b),  respectively.  The  yield  surface  is 
similar  to  that  calculated  at  zero  pressure,  and  can  be 
well  described  using  the  von  Mises-type  ellipse.  The 
center  of  the  ellipse  is  shifted  from  the  ideal  value 
(-1,-1)  by  a  very  small  amount  a  in  the  third  quadrant 
due  to  the  inherent  asymmetry  in  tension  and  compression. 


The  lengths  of  the  major  and  minor  axes  are  tabulated  in 
Table  II  along  with  the  shift  in  the  center  of  the  ellipse 
for  the  yield  and  flow  surfaces.  The  fitting  parameters 
for  the  yield  stress  locus  give  the  length  of  the  minor  axis 
to  be  1.2  GPa,  which  compares  very  well  with  the  value 
of  -1.22  GPa  predicted  by  a  traditional  von  Mises 
ellipse  using  Oy.  The  fitting  result  also  gives  the  length 
for  the  major  axis  to  be  2.25  GPa,  which  is  slightly 
larger  than  the  value  of  -2.12  GPa  predicted  by  a 
traditional  von  Mises  ellipse.  Similarly,  the  flow  stress 
ellipse  has  a  minor  axis  of  1.66  GPa,  which  agrees  very 
well  with  the  value  of  -1.70  GPa  predicted  by  a 
traditional  von  Mises  ellipse  using  Oy.  The  length  of 
the  major  axis  for  the  flow  stress  ellipse  is  calculated  to 
be  3.35  GPa,  which  is  slightly  larger  than  the  value  of 
-2.95  GPa  predicted  by  a  traditional  von  Mises  ellipse. 

D.  3-D  Yield  Surface 

A  full  3-D  yield  surface  can  be  obtained  by  superpo¬ 
sitioning  all  the  yield  and  flow  surfaces  on  top  of  each 
other  along  the  Z  axis.  The  superpositioning  in  two 
dimensions  for  the  yield  and  flow  surfaces  is  shown  in 
Figures  9(a)  and  (b),  respectively.  It  can  be  seen  that  the 
ellipses  are  symmetric  along  the  Oi  =  02  axis.  This 
symmetry  is  inherently  predicted  by  the  von  Mises  yield 
criterion.  The  minor  axis  of  the  yield  stress  ellipse, 
however,  appears  to  be  fluctuating.  This  variation  in  the 
minor  axis  of  the  calculated  ellipses  may  be  attributed  to 
the  fluctuation  in  the  yield  stress  values  calculated  using 
the  0.2  pet  offset  in  the  stress-strain  curves. The  full 
surface  in  three  dimensions  for  the  yield  stresses  and 
flow  stresses  is  shown  in  Figures  10(a)  and  (b),  respec¬ 
tively.  The  ellipses  together  can  be  seen  to  form  a 
cylinder  that  is  symmetric  around  the  hydrostatic  axis. 

The  von  Mises-type  criterion  is  typically  applied  to 
polycrystalline  metals  in  which  the  deformation  mech¬ 
anism  is  dominated  by  dislocation-induced  plasticity. 
Nanocrystalline  metals  in  the  inverse  Hall-Petch  regime 
show  a  significant  contribution  to  plastic  deformation 
from  GB  sliding.  The  results  discussed  here  suggest  the 
von  Mises-type  criterion  be  applicable  to  predicting 
yield  behavior  at  ultrafine  grain  sizes  for  deformation  at 
ultrahigh  strain  rates.  Deformation  at  ultrahigh  strain 
rates  limits  the  GB  sliding  contribution,  rendering  the 
dislocations  to  be  the  dominant  carriers  of  plasticity. 
The  small  deviation  in  the  aspect  ratio  of  the  fitted 
ellipses  can  be  attributed  to  the  small  contribution  from 
GB  sliding  that  may  still  exist. 


V.  CONCLUSIONS 

Large-scale  MD  simulations  are  carried  out  to  under¬ 
stand  the  macroscopic  deformation  behavior  of  nano¬ 
crystalline  Cu  with  a  grain  size  of  6  nm  at  high  strain 
rates.  Three  aspects  of  deformation  behavior  are  stud¬ 
ied:  the  tension-compression  strength  asymmetry,  the 
biaxial  yield  surface,  and  the  3-D  yield  surface.  The  MD 
simulation  results  suggest  that  strength  asymmetry  is 
sensitive  to  the  strain  rates  used.  Higher  strain  rates  give 
a  greater  asymmetry,  with  the  material  being  stronger  in 
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Fig.  9 — Superposition  of  {a)  yield  stress  (0.2  pet  offset  stress)  surface 
and  {b)  flow  stress  (peak  value  of  stress)  surface  at  various  values  of 
constant  Z  stress  (dj.  Ellipses  are  symmetric  along  the  oi  =  02  axis, 
as  predicted  by  the  von  Mises  yield  criterion. 

tension  as  compared  to  tension.  The  calculated  biaxial 
yield  surface  for  the  yield  and  flow  stresses  can  be  fit  to 
an  ellipse,  suggesting  the  applicability  of  a  von  Mises- 
type  yield  criterion.  The  aspect  ratio  of  the  calculated 
ellipse  differs  very  slightly  from  that  for  the  ellipse 
predicted  using  the  traditional  von  Mises  yield  criterion. 
These  results  indicate  that  the  von  Mises-type  yield 
criterion  may  be  used  to  study  deformation  at  high 
strain  rates  for  nanocrystalline  Cu  at  grain  sizes  in  the 
inverse  Hall-Petch  regime.  The  applicability  of  the  von 
Mises  yield  criterion  at  ultrahigh  strain  rates  is  attrib¬ 
uted  to  the  limited  GB  sliding  contribution  to  plastic 
deformation  at  ultrahigh  strain  rates.  Plasticity  is 
therefore  controlled  by  dislocations  at  high  strain  rates. 


(b) 


Fig.  10 — Plot  of  3-D  surface:  {a)  yield  stress  (0.2  pet  offset  stress) 
surface  and  {b)  flow  stress  (peak  value  of  stress)  surface  at  various 
values  of  constant  Z  stress  (cj.  Ellipses  form  a  cylinder  in  three 
dimensions  and  indicate  the  validity  of  the  von  Mises  yield  criterion 
for  plastic-yield  behavior  of  nanocrystalline  Cu  at  high  strain  rates. 


The  results  discussed  here  do  not  suggest  the  von  Mises 
criteria  for  quasi-static  deformation  of  ultrafine-grain- 
size  nanocrystalline  metals.  Significant  work  is  needed  to 
evaluate  the  contributions  of  GB  sliding  to  plasticity 
under  multiaxial  loading  conditions  and,  in  turn,  the 
yield  criterion  to  predict  the  macroscopic  behavior  of 
these  metals  at  quasi-static  loading  conditions. 
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